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Abstract 

It is shown that a supersymmetric SO (10) model extended with fermion singlets can accommodate 
Oh| the observed neutrino masses and mixings as well as generate the desired lepton asymmetry in 

concordance with the gravitino constraint. A necessary prediction of the model is near-TeV scale 
QJ \ doubly-charged Higgs scalars which should be detectable at the LHC. 

The observed neutrino masses constitute a compelling evidence of interactions beyond the Standard 
Model of particle physics and leave an impact in areas as diverse as astrophysics, cosmology, nuclear 

Q\ \ physics, and geophysics. The smallness of these masses finds a natural explanation in the see-saw 
mechanism [I], which requires a heavy Majorana (self-conjugate) neutrino. Such heavy neutrinos 

£T) • appear in grand unified theories (GUTs) based on SO(10), which incorporate quark-lepton unification 
and left-right symmetry [21 [3]. The wide disparity between the weak and unification scales in these 
models calls for a protection mechanism and supersymmetry (SUSY) is widely considered to be an 
attractive candidate. Further, such a model with a low SUSY scale leads to a unification of gauge 
couplings at high energies. These positive features have encouraged many explorations of the SUSY 
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SO(10) model. 



Another open problem, also of much interest, is the origin of the observed baryon asymmetry of the 
universe. Originally it was expected that baryon number violation inherent in GUTs will lead to this 
small asymmetry when heavy gauge (and/or Higgs) bosons decay while they are out of equilibrium. 
This hope was belied however since any primordial GUT-origin asymmetry will be totally diluted in 
the inflationary epoch. This has provided impetus to look for lower energy avenues for generating this 
asymmetry. An oft-chosen route is to generate a lepton asymmetry through the C and CP-violating 
out-of-equilibrium decay of heavy Majorana neutrinos. This is later converted to a baryon asymmetry 
through anomalous (B + L) violation, which is implicit in the Standard Model [HE]. 

It is but natural to ask whether the heavy Majorana neutrino which drives the neutrino mass see-saw 
can also generate the lepton asymmetry through its decay. This would have been truly economical. 

Hindrances to this programme arise from several directions, (a) The observed light neutrino masses 
require the heavy neutrino, which is right-handed (RH), to have a mass ~ 10 13 GeV. This sets the 
scale, M R , for the SU(2) r x SU(2) l x U(1) (b _ l) -> SU(2) L xU(l) Y gauge symmetry breaking, (b) 



Within the SO(10) GUT framework, the intermediate symmetry breaking scales are fixed through the 
Renormalization Group (RG) equations which reflect the gauge couplings' evolution with energy. In 
the simplest SO(10) GUT it is well-known that M R turns out to be ~ 10 16 GeV. (c) In a SUSY context 
there is an additional constraint, namely, to ensure that there is no overabundance of gravitinos in 
the universe. To maintain consistency with this, it has been demonstrated El [8] that the lepton 
asymmetry must be generated through the decay of a heavy neutrino whose mass does not exceed 
~ 10 7 ~ 9 GeV in order to prevent a washout, whereas leptogenesis through the canonical Type-I see- 
saw mechanism sets the lower bound 4.5 x 10 9 GeV. These conflicting requirements have acted as 
obstacles to a successful implementation of this attractive possibility. 

In this letter we propose a remedy for these maladies confining ourselves to the SUSY SO(10) GUT. 
If sterile - i.e., SO(10) singlet - leptons are introduced, one for each generation EH ED 02] , then 
a novel way can be found to meet the demands outlined in the previous paragraph. 

The uncharged fermions in this model, per generation, are the following: a left-handed neutrino v, a 
right-handed neutrino, N, and a sterile neutrino, S. For the three generation neutral fermion system, 
the mass matrix on which we focus is: 



M„ = (v N c S) 1 



where mu, Mx, Mx, and \i are all 3x3 matrix blocks. 
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It is not unreasonable to expect that the mass matrix in eq. ([I]) will alleviate the tension, summarised 
earlier, between light neutrino masses and adequate low-scale thermal leptogenesis. As discussed 
below, the double see-saw structure for the light neutrino masses, arising from eq. (pQ), also decouples 
it to some extent from low-scale leptogenesis; Mx,Mx, and fx appear in different fashions in the 
expressions. Utilizing an extension of the Minimal Supersymmetric Standard Model (MSSM) by the 
addition of RH neutrinos and extra fermion singlets, which results in a neutrino mass matrix of the 
structure of eq. ([I]), Kang and Kim [12] have found solutions to both the above issues. There, mo 
has been identified, as is done in the MSSM, with the charged lepton mass matrix. On the other 
hand, in the SO(10) model which is espoused here, quark-lepton symmetry [2] identifies the neutrino 
Dirac mass matrix rri£> with the up-quark mass matrix whose 33 element is nearly 100 times heavier. 
This, along with other GUT constraints, pose additional hurdles in addressing the problems in SUSY 
SO(10). 

We work in a basis in which the down-quark and charged lepton mass matrices are diagonal. This 
ensures that the entire mixings in the quark and lepton sectors can be ascribed to the mass matrices of 
the up-type quarks and the neutrinos, respectively. Using quark-lepton unification, the quark masses, 
and the Cabibbo-Kobayashi-Maskawa mixing angles, one therefore obtains ma, upto 0(1) effects due 
to RG evolution. 

We utilize spontaneous symmetry breaking of SUSY SO(10) with the Higgs representations 210,54, 
126 © 126, 16 16, and 10. By using the mechanism of D-Parity breaking near the GUT scale 
[13], the RH-triplet pair A R (1, 3, 1, -2) 0^(1,3,1,2), and the RH-doublet pair xr(1, 2, 1, -1) © 
Xr(1, 2, 1, 1) are treated to have masses at much lower scales compared to their left-handed counter- 
parts. Mx = Fxr, in eq. ([T]), is generated via the vacuum expectation value (vev) < xjj(1, 2, 1, —1) >= 
< Xr(1, 2, 1, 1) >= xr, where we take F to be a matrix with entries O(0. 1). 

Although we do not assign any direct vev to the RH-triplets, through a Ar exchange involving a 
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trilinear coupling in the superpotential, XArXrXr, an effective vev < A#(l, 3, 1, —2) >= vr = A — R 
is generated, resulting in the mass term Mx ~ / < Ar(1,3, 1,— 2) >, where / is a typical Yukawa 
coupling of Majorana type. If m^ R is around 1 TeV, which can be arranged by a tuning of the D-parity 
breaking term in the Lagrangian, the entries of Mx are O(10 n ) GeV. Without any loss of generality, 
Mx can be chosen to be diagonal. 

Notice that the SU(2) R x U(1) B -l symmetry breaks at the scale < A R (1, 3, 1, -2) > ~ 10 11 GeV 
while xr ~ 10 7 GeV. The states A^ and Re(A° R ) are eaten up as Goldstone bosons by the W R and 
Wft fields and A^ + and Im(A R ) survive as physical states with mass ~ 1 TeV. The Type-II seesaw 
contribution to the light neutrino mass matrix is damped out in this case because of the large masses 
of the left-handed Higgs triplet leading to m/j ~ 10~ 5 eV - 10 -6 eV |144 115] . Further, the vev of \l 
is zero or negligible. 

Block diagonalization of the mass matrix in eq. (pTJ) in the limit in which we are working {i.e., 
M N > M x » /i » m D ) leads to: 

M 2 M 2 
m v ~ -m D [M x - l ix{M x T )- 1 ] m D T , M s ~ fi - j^, M ~ M N + (2) 

where m„, Ms, and M are 3x3 matrices. The light neutrino masses are in a double see-saw pattern and 
// is determined once Mx is fixed. It may be noted that the mass matrix structure in eq. (p]) ensures 
that the type I see-saw contribution is absent and Mx remains unconstrained by the light neutrino 
masses. This freedom in Mx - a hallmark of the model - is vital to ensure adequate leptogenesis. 
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Figure 1: The tree and one-loop contributions to the decay of T\ that generate the lepton asymmetry. 



The eigenstates of Ms, which we denote by 2$, (i = 1, 2, 3), are superpositions of the sterile neutrinos 
S (predominant) and the right-handed ones N. These states are found to lie well-below 10 9 GeV, 
consistent with the gravitino constraint. In fact we show that this model allows successful leptogenesis 
at a temperature T ~ Mt — 5 x 10 5 GeV, which is nearly 4 orders below the maximum allowed value. 
Further, the singlet fermions decay through their mixing with the iVj which is controlled by the ratio 
Mx/Mn. The latter, which have masses O(10 n ) GeV and are off-shell, decay to a final Icj) state, 
where I is a lepton doublet and <f> the up-type MSSM Higgs doublet. This two-step process - for which 
tree and loop diagrams are depicted in Fig. [Tj (SUSY contributions are small) - results in a lepton 
asymmetry of the correct order. Because of the large value of Mx 3> Mx, a small S — N mixing 
results naturally in the Tj which in turn guarantees the out-of-equilibrium condition to be realized 
near temperatures T ~ Mt- 

A quantitative analysis of this programme has been carried out using the Boltzmann equations deter- 
mining the number densities in a co-moving volume Yt = nx/ns and Yl = n^jns, where ut, til an d 
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rig are respectively the number densities of the decaying neutrinos, leptons and the entropy: 
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§^ = eT ^rh ( Yt ~ y t) ~ -&rr^- ( 3 ) 

dz zH(z) ' zH(z) 



where T^, and T^y represent the decay, scattering, and wash-out rates, respectively, that take part 
in establishing a net lepton asymmetry. We refrain from presenting their detailed expressions here 
|16j . The Hubble expansion rate H(z), where z = My/T, and the CP-violation parameter are given 
by 

H{M T ) 1/2 M| T(T -+ lct>) - T(T -> fe*) 

^ (Z) = ^^' ff W = l-67^ — , ZT = T ^— WTW ^ Wy (4) 

Our target is to use eqs. ([2]) and (J3j) to obtain an acceptable solution within the framework of SUSY 
SO(10). Through an exhaustive analysis we find an appropriate choice of the block matrices appearing 
in eq. ([T]) which guarantees adequate leptogenesis while maintaining full consistency with the observed 
neutrino masses and mixing as well as the gravitino constraint. The mass scales are fixed as dictated 
by the RG evolution of gauge couplings in SUSY SO(10) when effects of two dim. 5 operators scaled 
by the Planck mass are included jl6L I17|, [T8] . The strategy we follow is to choose the matrix Mx first. 
To minimize the number of independent parameters, we take the matrix F to be real and diagonal, 
which is reflected in Mx- Then using mo, as fixed by quark-lepton unification, [i is determined from 
the double see-saw formula given in eq. ([2]). Using these inputs, one has to examine, by trial and 
error, different choices of Mn for adequate lepton asymmetry generation. 




Figure 2: The comoving density of T\ - Yt - and the leptonic asymmetry - Yl - as a function of z. Also shown 
is Y^ q . The inset displays the decay (r^ ) and inverse-decay (T T D ) rates of T\ compared with the Hubble expansion 
rate, H , as a function of z. 



The results for the development of the leptonic asymmetry as the universe evolves are shown in Fig. [2j 
They are obtained with the choice Mx = diag (0.2, 0.3, 0.4)xr with xr = 6 x 10 6 GeV. The neutrino 
Dirac mass matrix, mc, is constructed utilizing quark-lepton symmetry; the up- type quark mass 
eigenvalues and the CKM mixings are taken at the PDG [19] values with the CKM-phase as 1 radian. 
The neutrino masses are fixed so as to satisfy Am^ = 8.0 x 10 -5 eV 2 and Am| 2 = 2.5 x 10 -3 eV 2 
with the lightest neutrino taken massless. The neutrino mixing angles used are 623 = 45°, #12 = 32°, 
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and #13 = 7°. No other CP-phases are introduced in the lepton sector except the one through the 
CKM matrix for quarks. With these inputs the matrix /i is calculated following eq. ([2]). For the 
RH-neutrino we use the mass matrix = diag(0.1, 0.5, 0.9) x 10 11 GeV. This is consistent with 
m A R ~ 1 TeV. We assume that in the very initial stages the number densities, Y^, i = 1,2,3, and 
the leptonic asymmetry, Yj,, are zero. The chosen input values of the mass parameters result in a Tj 
mass spectrum such that only one state - T\ — is above the kinematic threshold for l<p production 
(rriTx = 3.9 x 10 5 GeV) and the lepton asymmetry results through its decay. This ensures that the 
leptogenesis is consistent with the gravitino bound. It is seen from Fig. [2] that T\ decays fall out of 
equilibrium as the universe expands (inset) and Yl achieves the right order (~ 10~ 10 ) starting off from 
a vanishing initial value while Yt steadily tends towards Y^ 9 . 

We stress again that an important outcome of the symmetry breaking is that out of the triplet 
Ar the components A^ and ReA R are absorbed as longitudinal modes of the broken generators of 
SU(2)r x U(1)b-l- The physical states are A R + and ImA° R and their superpartners. They will be 
within striking range of the LHC and the ILC with itia — 300 GeV - 1 TeV. 

Finally, we briefly discuss the mechanism of SUSY SO(10) breaking |X6 |, 117 } [IB]: 

50(10) SU(S) C x SU(2) R x SU(2) L x U(1) B -l [63221] 

^ SU(3) C x SU(2) L x U(1) Y [Gstd}^SU(3)cxU(l) Q (5) 

The SO(10) Higgs multiplets 210 and 54 are utilized to break the symmetry at Mjj. Within the 
210 there are two components which develop vevs; one breaks SO(10) to 63221 while the other is 
responsible for D-parity breaking. The vev of the singlet under the Pati-Salam group contained in 54 
ensures that there are no light pseudo-goldstone bosons arising from the 210 to upset perturbative 
gauge coupling evolution. As already discussed, SU(2)r x U(1)b-l is broken by the vevs of RH- 
triplets in 126 126. This induced vev, vr ~ 10 11 GeV, is also responsible for the masses of the iVj. 
The last step of breaking in eq. © relies on the electroweak vev of the weak bi-doublet in 10. We 
have carried out an analysis of the RG evolution of the gauge couplings to determine the intermediate 
mass scales. We find that Mr ~ 10 9-11 GeV can be obtained through the introduction of effective 
dim. 5 operators scaled by the Planck mass, Mpi [18]. It is noteworthy that both 210 and 54 are 
necessary for a viable SUSY SO(10) breaking pattern and that the resulting two dim. 5 operators are 
instrumental in alleviating the problem of leptogenesis under the gravitino constraint: 



The details of this analysis will be presented elsewhere |16j . Suffice it to state that \r]\^\ ~ 0(1) 
and the interactions in eq. (|6|) lead to finite corrections to the gauge couplings at the GUT-scale. 
The couplings of the left-right gauge group thus emerge from one effective GUT-gauge coupling. The 
upshot of this is that with these additional contributions it is possible to lower Mr to as low as 
10 9 — 10 11 GeV as required in this model. The grand unification scale is high: Mjj ~ 10 17 ~ 18 GeV 
and the model predicts a stable proton for all practical purposes. 

We expect that this model will have a natural extension to an E(6)-GUT wherein the matter multiplets 
and the singlet fields will constitute the fundamental 27 representation of the gauge group. 

In conclusion, we have presented a SUSY SO(10)-based model relying on a double see-saw mechanism 
which is (a) consistent with the known neutrino masses and mixing, and (b) can lead to a correct lepton 
asymmetry via the decays of sterile, i.e., SO (10) singlet, neutrinos while remaining in concordance 
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with the gravitino constraint. The intermediate scales are obtained through an RG analysis of the 
gauge coupling running and are consistent with a long-lived proton. The model is falsifiable through 
its prediction of doubly-charged Higgs bosons within the reach of the LHC. 
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